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PREFATORY  NOTE 


The  casual  reader  may  not  be  familiar  with  the  term  "waste 
heat"  as  used  here  or  with  the  qualitative  aspects  of  the  temperature 
history  of  a  material  element  which  is  shocked  and  then  allowed  to 
expand  to  ambient  pressure.  He  can  find  these  ideas  discussed  in 
Chapter  IX  of  Physics  and  Chemistry  of  High  Pressure,  edited  by 
R.  S.  Bradley  and  published  by  Academic  Press,  196?,  and  in  some 
of  the  references  given  there  to  other  papers  by  G.  E.  Duvall.  If  he 
does  not  wish  to  go  so  far,  the  following  description  may  be  helpful. 

The  process  of  shock  compression  is  necessarily  dissipative, 
so  the  temperature  of  a  material  element  elevated  to  a  high  pressure 
state  by  passage  of  a  shock  is  necessarily  higher  than  the  initial  tem¬ 
perature  of  the  unshocked  material.  Call  this  elevated  temperature 
the  "shock  temperature"  of  the  material  corresponding  to  the  shock 
pressure  P  .  This  increase  in  temperature  results  from  the  dissi¬ 
pative  forces  coming  into  play  in  the  shock  process  and  also  from 
isentropic  compression,  aside  from  any  dissipation.  The  material 
element  is  returned  to  ambient  pressure  (whence  it  started  before 
shock  compression)  by  a  rarefaction  process  in  which  dissipative 
forces  can  be  neglected.  The  isentropic  expansion  cools  the  material 
to  a  final  temperature  which  is  lower  than  the  shock  temperature  but 
higher  than  the  initial  temperature  of  the  unshocked  state.  Call  this 
the  "final  temperature"  corresponding  to  the  ambient  pressure  PQ  . 
Since  the  final  temperature  is  greater  than  the  initial  temperature, 
and  the  pressures  for  initial  and  final  states  are  the  same,  the  passage 
of  the  shock  and  subsequent  rarefaction  have  increased  the  specific  en¬ 
thalpy  of  the.  material  element.  This  increase  in  enthalpy  is  called 
here  the  "waste  heat"  deposited  by  the  shock  in  unit  mass  of  the  mate- 


ABSTRACT 


Procedures  for  relating  waste  heat  to  entropy  increases  in  a 
shockwave  are  described.  Numerical  coefficients  for  Mie -Gruneisen 
equations  of  state  are  obtained  for  several  liquids,  assuming  that 
(BP/dT)v  and  Cv  are  constant.  These  equations  of  state  are  used 
to  calculate  the  flow  produced  by  a  1.  b-cm-radius  sphere  of  C-3-like 
explosive  in  a  surrounding  liquid.  Specific  entropy  increase  as  a 
function  of  radius  is  included  in  the  calculation,  and  this  is  used  to 
determine  the  waste  heat  for  a  3.  85 -cm-radius  liquid  burster.  (Liquid- 
to -explosive  mass  ratio  is  about  10. )  Calculations  for  glycerine,  for 
two  different  water  equations  of  state,  and  for  ether  yield  total  waste 
heats  equal  to  (.32  ±.03)  x  (explosive  energy).  Close  agreement 

between  calculations  for  two  different  water  equations  of  state  pro¬ 
duces  considerable  confidence  in  the  reliability  of  this  result.  On  the 
basis  of  these  calculations,  evaporation  of  several  liquids  is  estimated. 
If  decomposition  can  occur  or  if  more  precise  values  of  waste  heat  are 
needed,  more  exact  equations  of  state  will  be  required;  the  present 
calculations  can  be  readily  modified  to  accommodate  such  equations. 
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I  INTRODUCTION 


Flash  radiography  of  the  expansion  of  a  liquid-filled  spherical 
burster  and  calculations  of  the  kinetic  energy,  energy  of  drop  forma¬ 
tion,  and  energy  of  rupture  of  the  plastic  case  indicate  that  a  large 
fraction  of  the  explosive  energy --perhaps  50%--is  not  accounted  for 
by  these  mechanisms.  Part  of  this  is  dissipated  in  the  liquid  of  the 
fill  as  heat  because  of  the  inherent  dissipative  nature  of  shock  com¬ 
pression.  It  is  the  purpose  of  this  calculation  to  determine  what 
fraction  of  the  explosive  energy  is  dissipated  in  this  way  for  a  parti¬ 
cular  burster  and  to  outline  a  procedure  for  calculating  the  energy  so 
dissipated  in  a  more  general  case.  To  this  end  we  proceed  to  formu¬ 
late  equations  of  state  for  a  number  of  liquids  and  to  employ  a  spheri- 
cal  Q-code  to  calculate  flow  produced  in  the  liquid  by  detonation  of 
the  explosive  and  the  dissipated  energy. 


We  shall  call  this  dissipated  energy  the  "waste  heat."  An 
expression  for  it  is  required  which  can  be  used  in  the  flow  calculation, 
given  a  suitable  equation  of  state.  In  the  shock  process  an  element  of 
mass  is  compressed  to  a  state  (Pj,  Vj)  and  then  it  expands  along  an 
isentrope  to  the  initial  pressure  PQ  and  to  a  volume  ,  greater 
than  the  initial  volume  VQ  (see  Fig.  1).  The  change  in  enthalpy  in 
going  from  the  initial  state  (PQ,  V  )lo  the  final  state  (PQ ,  V^)is 


simply 


See,  for  example,  "Computational  Methods  in  Physics,  "  edited  by 
Fernbach  and  Alder,  Vol.  Ill,  "Hydrodynamics,"  Academic  Press, 


FIG.  1  P-V  DIAGRAM 


where  Tq  and  Sq  are  temperature  and  entropy  of  the  initial  state 
(PQ,  VQ)  and  (T^,  )  correspond  to  (PQ,  V^).  The  change  in 

internal  energy  corresponding  to  the  expansion  is 

E  '  -  E  *(H#-H)-P(V'-V);  (2) 

o  o  '  o  o'  o  o  o  ' 

H'  -  HL  is  the  waste  heat  per  unit  mass.  Since  P  ~  0  ,  H  *  -  H  s 
o  o  r  o  o  o 

E^  -  E^  to  an  excellent  approximation.  In  order  to  determine  E^ , 

,  or  directly,  the  flow  calculation  must  be  continued  until  the 
pressure  returns  to  PQ  ,  or  very  close  to  it.  This  is  an  impractical 
procedure  with  a  Q-code  for  two  reasons:  first,  the  computing  time 
would  be  very  long,  hence  expensive.  Second,  the  round-off  errors 
creeping  into  the  calculation  would  be  apt  to  give  serious  errors  in 
the  end  result. 

An  alternative  and  more  practical  procedure  is  to  calctilate  the 
entropy  produced  by  passage  of  the  shock.  Since  it  increases  during 
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shock  compression  and  remains  constant  in  the  subsequent  rarefaction, 
it  can  be  calculated  during  that  part  of  the  numerical  integration  where 
the  precision  is  high;  when  it  is  known  everywhere  within  the  volume 
of  interest,  the  waste  heat  can  be  calculated  directly.  We  have  on  the 
isobar  P  =  PQ  , 

C  dT 

dS  =  ,  (3) 


and  the  indefinite  integral  gives  S(T).  If  this  can  be  solved  for  T(S), 
the  solution  can  be  substituted  into  the  second  integral  in  Eq.  (1)  and 
is  then  given  explicitly  in  terms  of  S^.  If  Eq.  (3)  cannot  be 
solved  for  T(S),  it  is  necessary  to  resort  to  a  numerical  procedure 
in  which  the  second  term  on  the  right-hand  side  of  Eq.  (E)  is  evaluated 

from  the  integral  , 

o 

Vo-  Vo  =  /  OpdT-  <4> 

To 

using  known  values  of  the  thermal  expansion  coefficient. 


For  the  particular  case  C  =  constant  =  C  and  («)  =  constant, 

p  v  oT  p 


we  have: 


So-So  =  VnTo/To 
To  =  ToexPt<So-So>/Cp3 


(5) 


Vo-  V0  =  <IS)«(T>TJ 


Fpuo  O' 

(  ItV  To  {expt(So‘ So)/CpJ  " 


(6) 
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-  Eq  =  Cp(T^  -  T0)  -  PQ  (V^  -  Vc) 

(7) 

=  tcP-po<lr»p]To{exp[(so'V/cp]  •  *} 

For  the  liquids  listed  in  Table  I,  («)  <  3  x  10~4cc/g°C  and 

~6  o  P  _ ^ 

C  >  8  x  10“  mbar  cc/g°C.  Then  since  P  =1  atmosphere  =  10 
p  o 

mbar,  we  have 

<  10”4  (8) 

P 

Consequently,  to  a  good  approximation  -  Eq  =  -  Hq  and,  in  this 

case, 

Eo  -  Eo  *  CpTo  {exp  1  <So  -  So>/Cp]  -  1  }  (9> 

Then  if  the  liquid  surrounding  the  explosive  is  partitioned  into  spheri¬ 
cal  shells  of  mass  Kij  ,  and  -  SQa  is  calculated  for  each  cell, 
we  have  for  the  total  waste  heat  in  J  cells  behind  the  shock: 

J 

W  =  CnTo  Z  m  [exp (S./C  )  -  l]  (10) 

r  W  J  J  Sr 

j  =  l 

We  shall  use  this  expression  for  calculating  waste  heat  for  the  present 

gy 

problem.  The  assumptions  that  Cp  =  Cv  =  constant  and  that  )p  ~ 
constant  are  unpalatable,  but  there  is  some  reason  to  suppose  that 
waste  heat  is  quite  insensitive  to  assumptions  about  specific  heat  any¬ 
way.  Consider  the  following  argument:  For  very  weak  shocks 
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II  EQUATIONS  OF  STATE 


A.  General 

An  equation  of  state  of  a  nonreactive  flud  is  a  relation  between 
three  thermodynamic  variables  which,  given  two,  allows  the  third  to 
be  calculated.  Such  an  equation  may  be  incomplete,  i.e.,  inadequate 
for  calculation  of  other  thermodynamic  variables.  For  example,  the 
ideal  gas  equation  of  state  is  a  relation  between  P,  V,  and  T,  but  the 
internal  energy  cannot  be  calculated  without  further  information.  On 
the  other  hand,  a  relation  between  E,  S,  and  V  is  complete  since 

<!i>v  =  T<s-  v> 

<H>„  -  -p<s-v> 

H  =  E  +  PV  =  H(S,  P) 

F  =  E  -  TS  =  F(T,  V) 

G  =  H  -  TS  =  G(T,  P) 

In  order  to  construct  an  equation  of  state  it  is  necessary  to  fit 

whatever  experimental  data  are  available  into  a  frame  which  is  thermo¬ 
dynamically  consistent  and  which  does  not  violate  our  concept  of  what 
a  "well-behaved"  material  should  do. 


These  remarks  are  suggested  by  conversations  with  Dr.  Michael  * 

Cowperthwaite  of  Poulter  Laboratories  and  by  his  memorandum: 
"Comments  on  the  Calculation  of  Shock  Temperature  along  Hugoniot 
Curves." 
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Available  data  may  coneiat  of  thermal  expansion  coefficients  vs. 
temperature  at  atmospheric  pressure,  sound  velocity  at  atmospheric 
pressure,  P-V  data  obtained  statically  under  isothermal  conditions, 
and/or  shock  compression  data.  Occasionally  data  on  expansion  co¬ 
efficients  at  elevated  temperatures  and  pressures  are  available.  It 
is  further  assumed  here  that  isotherms  and  adiabats  of  well-behaved 
materials  satisfy  the  condition  (g^)*j>orS  <  ®  everywhere  and 
(-~2  )xorS  >  ®  everywhere.  We  shall  ignore  phase  changes  in  what 

follows,  except  for  a  discussion  of  vaporization  in  Section  V. 

B.  Gruneisen  Equation  and  Specific  Heat 

the  calculation  of  wave  propagation,  a  (P,  V,  E)  equation 
of  state  adequate,  and  we  shall  assume  the  Mie -Gruneisen  form: 


P(V,  E)  =  P.(V)  +  £[E  -E.(V)]  , 


where  P^(V)  and  E^(V)  are  pressure  and  internal  energy  on  the  iso¬ 
therm  passing  through  (P  ,  V  ),  and  T  =  r(V)  is  the  Gruneisen  ratio 

1  ®  O  I 

r  =  aV/0Cv,  where  a  =  -^  (^j)p  end  3  =  -  y  •  Some  values  of 

T_  (V  ,  P  )  and  T  /V  are  shown  in  Table  I. 
o  o  o  o  o 

For  metals,  where  a  theory  of  the  variation  of  T  with  V 
exists,  *  it  is  found  that  r/V  =  constant  is  a  fair  approximation  to 
theoretical  behavior.  That  assumption  will  be  made  here  and  the 
ratio  denoted  by  the  symbol  b  .  Values  of  the  specific  heat  are  re¬ 
quired  in  order  that  entropy  be  calculable,  and  the  simplest  possible 
assumption  is  also  made:  Cy  =  =  constant.  This  condition  violates 

our  desire  for  thermodynamic  consistency.  The  ratio  C^/Cv  is  given 
by  the  expression 

Q 

P  _  l  i  m  T  /  3V  \  _  i  i  v  T>  / 1  \ 


1  +  tV(StV  1+bVaT  . 


Rice,  McQueen,  and  Walsh,  Solid  State  Physics,  Vol.  6,  pp  1,  ff, 
Seitz  and  Turnbull,  Eds.,  Academic  Press,  1958. 
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From  the  entries  in  Table  I  it  can  be  found  that  b  VaT  is  as  large  as 
0.  6  at  20° C,  1  atmosphere,  so  that  the  assumption  Cy  =  Cp  is  not 
a  particularly  good  one.  However,  variations  in  Cp  as  well  as  un¬ 
certainty  about  its  average  value  over  an  extended  range  of  pressures 
and  temperatures  are  such  that  this  assumption  is  in  keeping  with 
other  uncertainties  in  the  equation  of  state.  We  thus  arrive  at  the 
following  equation  of  state: 

P  =  P.(V)  +  b[E  -  E.(V)]  (13) 

E  -E.(V)  =  Cv(T-  Tq)  .  (14) 

The  amount  of  time  required  to  formulate  an  equation  of  state 
which  is  much  better  than  that  represented  by  Eqs.  (13)  and  (14)  would 
be  substantial,  and  it  is  thought  that  the  increase  in  accuracy  thereby 
obtained  does  not  now  warrant  the  additional  investment.  If  the  results 
of  the  present  calculation,  concerning  waste  heat,  temperature  rise, 
and  vaporization,  point  to  possible  serious  consequences  in  burster 
performance,  a  more  careful  study  may  be  justified. 

C.  Evaluation  of  P^(V) 

Since  we  are  concerned  with  pressures  of  several  kilobars  or 

more,  we  can  set  PQ  =  0  without  affecting  results,  and  we  do  so 

here.  The  data  to  be  used  in  determining  P^(V)  are- the  values  of 

b  and  C  in  Table  I,  points  on  the  Hugoniot  (R-H)  curve  measured  by 
v  2 

Walsh  and  Rice,  (cf.  Table  II),  and  values  of  adiabatic  compressibility 
obtained  from  sound  velocities  at  atmospheric  pressure.  For  points 


John  M.  Walsh  and  Melvin  H.  Rice,  J.  Chem.  Phys.  26,  815-823 
(1957). 
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Table  II 


EXPERIMENTAL  HUGONIOT  CURVE  DATA  FOR  LIQUIDS* 
(Velocities  in  km/sec,  pressure  in  kilobars,  VQ  in  cc/g(water 
data,  for  which  actual  initial  temperatures  were  15°C  to  30°C 
are  corrected  to  the  uniform  initial  state  T0  =  20°C,  PQ  =  1 
atmosphere)}. 


Liquid 

T 

(°C) 

V 

o 

U 

s 

U 

P 

P 

v/v 

o 

Water 

20 

1.0018 

3.354 

0.  952 

31. 

8 

0.716 

20 

1.0018 

4.093 

1.392 

56. 

8 

0.  660 

20 

1.0018 

4.  126 

1.411 

58. 

2 

0.  658 

20 

1.0018 

4.536 

1.655 

74. 

9 

0.  635 

20 

1.0018 

4.813 

1.829 

87. 

8 

0.  620 

20 

1.  0018 

4.777 

1.806 

86. 

1 

0.622 

20 

1.0018 

4.757 

1.798 

85. 

3 

0.  622 

20 

1.  0018 

5.626 

2.385 

133. 

9 

0.576 

20 

1.0018 

5.604 

2.370 

132. 

5 

0.577 

20 

1.0018 

5.601 

2.335 

130. 

5 

0.583 

20 

1.0018 

8.07 

4. 13 

330. 

0 

0.488 

20 

1. 0018 

8.  07 

4.24 

342. 

0 

0.  476 

20 

1.0018 

8.45 

4.60 

388. 

0 

0.  456 

20 

1.0018 

8.49 

4.72 

400. 

0 

0.444 

20 

1.0018 

8.59 

4.72 

405. 

0 

0.450 

20 

1.0018 

8.74 

4.81 

419. 

0 

0.450 

Mercury 

25 

0. 0739 

2.752 

0.608 

226. 

4 

0.779 

17 

0.  0739 

3. 101 

0.772 

324. 

0 

0.751 

24 

0.  0739 

3.504 

0.978 

463. 

7 

0.  721 

Acetone 

26 

1.274 

5.37 

2.510 

105. 

8 

0.533 

30 

1.279 

3.97 

1.495 

46. 

4 

0.  623 

Benzene 

32 

1. 155 

5.66 

2.470 

121. 

0 

0.564 

16 

1. 133 

4.  10 

1.448 

52. 

4 

0.647 

Bromoe  thane 

19 

0.  685 

4.  68 

2.300 

157. 

1 

0.508 

16 

0.  682 

3.40 

1.353 

68. 

0 

0.599 

Carbon 

33 

4.32 

2.412 

129. 

5 

0.441 

Disulfide 

17 

3.37 

1.415 

58. 

5 

0.580 

Carbon 

25 

0.  634 

4. 85 

2.235 

171. 

0 

0.539 

Tetrachloride 

22 

0.  629 

3.51 

1.325 

73. 

9 

0.  622 

Ethyl  Ether 

32 

1.433 

5.40 

2.550 

96. 

1 

0.  528 

H 

1.407 

3.88 

1.517 

41. 

8 

0.  609 

Table  II  (Concluded) 


Liquid 

T 

o 

(°C) 

Vo 

U 

s 

U 

P 

P 

V/V 

o 

Ethyl 

26 

1.275 

5.  63 

2.500 

110.4 

0.556 

Alcohol 

21 

1.267 

4.03 

1.487 

47.3 

0.  631 

Glycerine 

30 

0.798 

6.07 

2.240 

170.3 

0.631 

18 

0.794 

4.  58 

1.328 

76.6 

0.710 

Hexane 

32 

1.499 

5.54 

2.590 

95.7 

0.533 

19 

1.471 

4.02 

1.517 

41.5 

0.622 

Methanol 

24 

1.271 

5.51 

2.525 

109.5 

0.  542 

15 

1.255 

3.95 

1.483 

46.6 

0.  625 

N-amyl 

23 

1.236 

5.81 

2.465 

115.9 

0.576 

Alcohol 

19 

1.227 

4.26 

1.466 

50.9 

0.  656 

Toluene 

19 

1.  138 

5.73 

2.412 

121.5 

0.579 

El 

1.  141 

4.  12 

1.443 

52.  1 

0.650 

Mononitro- 

0.856 

5.  64 

2.300 

151.5 

0.592 

toluene 

H 

0.856 

4.20 

1.340 

65.8 

0.  681 

3$e  1 

This  is  Table  I  of  Reference  2. 
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P  (V)  on  the  R-H  curve,  Eq.  (13)  becomes 
H 


where 


P«  =  Pi+b(EH.E.) 


pu<V„  -  V)  +  E 


'H  "  2T  H '  o 


05) 

(16) 


Ei  =  Eo  +  / 


(17) 


=  E  +  b  C  T  (V  -  V  ) 
o  VO  o 


T 

-  /  pidv 


(18) 


Elimination  of  E„  and  E.  from  Eqs.  (15),  (16),  and  (18)  yields 
ri  l 

V 

PH  =  Pi  +  lPH  (Vo  *  V)  -  b2cvTo  ( V  -  VQ)  +  b  J  PjdV 


(19) 


Differentiation  of  Eq.  (19)  yields  the  following  differential  equation  for 


P.  in  terms  of  P„  and  its  derivatives: 
l  H 


dPi  .  ,  „  r,  b,„  dPH  ,  b~  .  .2, 


~27 


+  bpi  =  [  1  -  2-(Vq  -  V)]  -gy-  +  ?PH  +  b  CvTo 


(20) 


The  required  solution  of  Eq.  (20)  is 


P.[(V)  =  exp(- bV)]  [f(V)  -b  J  f(V)dV]  +bCvTo{l  -  exp[b(VQ- V)]} 


(21) 


where 


f(V)  =  l  1  -  \  (VQ  -  V)]  PH  exp  (bV) 


(22) 
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The  integral  in  Eq.  (21)  is  readily  evaluated  if  f(V)  is  fitted  by  a  poly¬ 
nomial  in  V  .  A  suitable  expression  is 

i{V)  =  f(Z)  =  AZ+A1Z2+A2Z3  (23) 

where 

Z  =  1  -  V/VQ  .  (24) 

Substitution  of  Eq.  (23)  into  Eq.  (21)  yields  the  following  expression 
for  P.(V): 

P.(V)  =  exp(-bV)  (AZ  +  BZ2  +  CZ3 +DZ4)  +  bC  T  [l-exp(bV  Z)] 

i  v  °  °  (25) 

where 

B  =  A,  +  bV  A/2 
1  o 

C  =  A2+bVoAj/3  (26) 

D  =  bV  A,/4 
o  c 


Since  f(V)  is  defined  along  the  R-H  curve,  its  slope  at  atmos¬ 
pheric  pressure  corresponds  to  adiabatic  compressibility  obtained 
from  the  sound  velocity  cq  .  This  determines  the  coefficient  A  : 


(—  ) 

'  dV 'V  =  V 

o 


A 


exp(bV0)  (-37  )y  = 
o 

C2  exp(bVo)/Vo  . 


co2exp(bVo) 


(27) 


Values  of  A  for  several  liquids  are  entered  in  Table  III. 


Table  III 


R-H  DATA  CORRECTED  TO  20°C 
AND  COEFFICIENTS  FOR  EQ.  (25) 


Liquid 

v/v 

o 

_ 

PH 

(Mbar) 

A 

B 

C 

D 

(Mbar) 

Acetone 

m 

.  1051 
.  0458 

. 03036 

-.2465 

1.  676 

.4301 

Benzene 

.  6470 
.5640 

.  0527 
.  1191 

. 05336 

-.2713 

2.412 

.  7991 

Carbon  Tetra¬ 
chloride 

.  6220 
.  5390 

.  0737 
.  1699 

. 05415 

-.4949 

3.  121 

1  .  110 

Ethyl 

Alcohol 

.  6310 
.  5560 

.  0473 
.  1100 

.  02008 

-.  7433 

2.965 

.4943 

Ethyl 

Ether 

.  6090 
.5280 

.0418 
.  0951 

. 01510 

-.4415 

1.872 

.  3670 

Glycerine 

.  7100 
.  6310 

.  0766 
.  1703 

.09914 

-.4080 

5.  088 

1.  010 

Methanol 

.  6250 
.  5420 

.  0468 
.  1093 

. 01812 

-.5231 

2.  320 

.  3641 

Water 

.  7160 
.  5760 

.  0318 
.  1329 

. 02443 

-.4433 

2.  723 

.07339 

Mercury 

.  7790 
.  7210 

.  2256 
.4626 

4.432 

.3767 

53.40 

40.22 

I 


Except  for  mercury  and  water,  there  are  but  two  points  on  the 
R-H  curve  for  each  liquid.  These  junt  suffice  to  determine  the  co¬ 
efficients  Aj  and  A ^  of  Eq.  (23).  Before  this  is  done,  the  data  should 
be  corrected  to  a  common  initial  temperature,  since  the  actual  temper¬ 
atures  vary  significantly  from  shot  to  shot.  This  can  be  accomplished 
by  finding  the  derivative  of 

PH(V’  Vo>  *  Pi<' v-  Vor>  +  I  PH  <' V  V)  +  bEo(To>  •  bEi(V>  V  •  (28> 


Here  the  isotherm  is  supposed  to  be  a  fixed  curve  passing  through  the 

reference  volume  V  on  which  T  =  T  .  Then  only  P.,,  E  .  and  V 

or  or  7  H  o  o 

are  functions  of  Tq  .  Differentiating  Eq.  (28)  gives  the  expression: 


dP._r  .  bPj.  dV  dE  bP„ 

<yr  M 1  -  l(vo-v»  ■  -r  ar +  b  dr2  ■  x  “.V  bc, 

o  o  o 


(29) 


Then  the  required  correction  to  P„  ,  in  order  to  refer  it  to  T  „ 

n  or 

ap. 

/  m  m  \  _  O  tt 

T 
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k  p«aV  +2C 

APH  =  (  jf'v  (Tor  ■  To*  =  2  .  b  vf  <Tor 

o  1  -  w(V  -  V) 


To»- 


(30) 


P..(V,  Tor)  was  calculated  for  the  liquids  listed  in  Table  HI.  Two  points 
each  were  selected  for  water  and  mercury  and  they  were  treated  as  the 
others.  Water  was  chosen  because  very  good  equations  of  state  exist 
for  it,  and  it  appeared  worthwhile  to  compare  the  simple  equation 
developed  here  with  a  more  elaborate  equation  of  state  developed  at 
U.S.  Naval  Ordnance  Laboratory.  Mercury  was  included  for  compari¬ 
son  with  other  metals.  Values  of  Ppj  corrected  to  20°C  are  given  in 
Table  HI. 


Coefficients  Aj  and  A^  of  Eq.  (23)  can  now  be  determined  by 
substituting  values  of  P^  and  V/VQ  from  Table  HI  into  Eq.  (23)  and 
solving  the  resulting  algebraic  equations  for  A^  and  A ^  •  These,  in 
turn,  substituted  into  Eqs.  (26),  yield  the  values  of  B,  C,  and  D 
entered  in  Table  HI. 
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Equations  (13)  and  (14),  the  data  entered  in  Table  II,  Eq.  (25), 
and  the  coefficients  in  Table  HI  represent  the  required  equation  of 
state  for  the  liquids  listed  in  Table  III, 
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HI  COMPUTED  ISOTHERMAL  AND  HUGONIOT  CURVES 


Values  of  pressure  (I3-)  and  energy  (E^)  on  the  293°K  isotherm 

and  values  of  pressure  (P,.),  energy  (E„),  temperature  (T„),  and 

ri  ri  ri 

entropy  (S^  -  So)  on  the  R-H  curve  passing  through  P  =  0,  T  =  293  K 
have  been  computed  from  Eqs.  (13),  (14),  and  (25)  with  the  aid  of  the 
Burroughs  5500  computer.  Values  of  these  quantities  as  functions  of 
specific  volume  are  entered  in  Table  IV  for  the  materials  listed  in 
Table  HI. 

Negative  values  of  PT. ,  Err,  S,.  -  S  ,  and  P.  in  Table  IV  are 

ri  ri  ri  O  l 

physically  meaningless;  they  are  a  consequence  of  fitting  f(V)  in 
Eq.  (23)  with  a  cubic  in  V  .  Essentially  the  same  result  would  have 
been  obtained  if  P^  had  been  fit  by  a  cubic.  The  inflection  point  occurs 
between  the  zero-pressure  point  and  the  first  shock  point,  and,  in  this 
region,  Eq.  (25)  is  not  a  good  representation  of  the  isotherm.  This 
misrepresentation  will  produce  some  error  in  energy,  entropy,  and 
temperature  parameters  calculated  at  higher  pressures,  but,  except 
for  ethyl  alcohol  and  ethyl  ether,  the  error  is  not  apt  to  be  serious. 

Negative  values  of  Ei  in  Table  IV  may  be  incorrect  in  magnitude, 
but  they  have  the  right  sign.  This  follows  from  Eq.  (18),  which  is 
differentiated  to  yield  the  equation 

dE.  =  (-  dV)  (P.  -  b  C  T  ) 

Here  (~dV)  >  0  ,  and  until  P^  exceeds  bCvTQ  ,  dE.  remains  negative. 

On  comparing  Ppj  for  water  with  the  Walsh  and  Rice  Hugoniot 
points  (Table  H),  it  is  apparent  that  between  31  and  130  kbars  the  com¬ 
puted  values  of  P^  are  too  high  but  do  not  differ  from  the  measured 
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points  by  more  than  10%.  Comparison  of  calculated  entropies  and 

3 

temperatures  with  those  given  by  Rice  and  Walsh  shows  the  tempera¬ 
ture  increases  to  be  half  the  Rice  and  Walsh  values,  while  the  entropy 
increases  are  about  half  again  as  great.  These  differences  are  princi¬ 
pally  due  to  the  high  value  of  specific  heat  used  here  compared  to  that 
used  by  Walsh  and  Rice  (1  cal/gdeg  compared  to  0.  86  at  25  kbar)  and 

the  variation  in  C  /(3V/3T)  allowed. 

P  P 


M.  H.  Rice  and  Johr  M.  Walsh,  J.  Chem.  Phys.  26,  824-830  (1957). 
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IV  WASTE  HEAT  CALCULATION 


The  flow  produced  in  an  unbounded  liquid  surrounding  a  spherical 
explosive  charge  has  been  calculated  numerically  by  direct  integration 
of  the  flow  equations,  using  a  Richtmyer-von  Neumann  artificial  vis¬ 
cosity  to  integrate  through  the  shock  front.  An  abbreviated  flow  chart 
showing  the  principal  steps  in  the  computation  is  given  in  Fig.  2. 

Some  details  of  the  computation  are  discussed  in  the  following  para¬ 
graphs.  The  code  used  in  the  calculation  was  developed  by  John  Erkman 
of  Poulter  Laboratories. 


Geometry 


The  radius  of  the  explosive  charge  was  taken  to  be  1.  6  cm. 
Detonation  was  initiated  at  the  center.  The  outer  radius  of  the  liquid 
spherical  shell  was  taken  to  be  4,0  cm.  The  computation  was  stopped 
before  the  wave  front  reached  this  outer  boundary,  so  the  liquid  was 
in  effect  boundless. 


Explosive  Properties 


Initial  density,  energy  of  explosion,  and  detonation  velocity  were 
taken  to  be,  respectively,  1.59  g/cc,  0.  05136  Mbar  cc/g,  and  0.8163 
cm/psec;  explosive  gases  were  assumed  to  be  polytropic  with  Y  =  2.736. 
This  loading  corresponds  approximately  to  a  low -density  RDX  loading; 
it  represents  a  total  energy  of  33,  500  calories,  compared  to  32,300 
for  Composition  B.  The  difference  is  small  and  in  the  direction  to  be 
expected  when  C-3  replaces  Composition  B. 


mation  of  State 


The  form  of  the  equation  of  state  and  the  parameters  were  dis¬ 
cussed  in  Sections  II  and  III.  The  computations  were  performed  for 
water  and  for  glycerine.  The  flow  chart  for  the  subprogram  which 
calculates  E  and  P  is  shown  in  Fig.  3. 
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FIG.  2  FLOW  CHART  FOR  WASTE  HEAT  CALCULATION 
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In  order  to  check  the  sensitivity  of  the  calculation  to  the 

assumptions  made  in  the  equation  of  state,  a  second  calculation  was 

made  for  water,  using  the  NOL  equation  of  state  as  programmed  by 

John  O.  Erkman.  It  is  believed  to  be  an  accurate  P-V-E  equation  of 

state  and  was  given  to  SRI  by  H.  Sternberg  of  the  U.S.  Naval  Ordnance 

Laboratory  in  Silver  Spring.  This  equation  was  supplemented  with 

specific  heat  and  thermal  expansion  coefficients  taken  from  Rice  and 
3 

Walsh.  It  probably  comes  close  to  being  the  best  equation  of  state 
for  water  that  can  be  constructed  with  presently  available  information. 
The  computing  procedure  is  quite  long,  but  it  is  reproduced  in  Fig.  4 
for  reference. 

D.  Waste  Heat 

This  calculation  is  made  in  the  main  program  after  the  shock 
front  has  advanced  to  a  radius  equal  to  the  outer  radius  of  the  liquid 
fill:  3.85  cm  in  the  present  case.  At  this  time  the  shock  calculation 
is  stopped,  the  waste  heat  is  calculated  for  each  cell  and  summed  for 
all  cells.  The  flow  chart  for  the  calculation  is  shown  in  Fig.  5. 
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FIG.  4  COMPUTATION  OF  E,  P,  T,  S,  USING  NOL  EQUATION  OF  STATE  OF  WATER 
AND  VALUES  OF  (<?V/dT)p,  Cp  FROM  RICE  AND  WALSH  (ibid.) 


5.  722427  x  10'3  -  1.240522Ej  +  50.42535  E^ 


-  1.400579  x  103E3  +  4.  13795  x  106e“J  -  2.  726437  x 

-  1.295684  x  101  1  E j  +  1.437988  x  1013Ej 

dFj/dEj 


FIG.  4  Continued 


j  Fz  =  -  2.74818  x  10"2  +  1.69113EJ  ♦  17.  12981E2 

♦  1.483384  x  104E3  -  1.549072  x  107E4 
I  1  l 

+  3.415519  x  109E*  -  2.357818  x  10U  E * 
F22’  dF2/dE» 


F2  =  -  2.21543  x  10"2  +  1.  51099E  -  10. 56399  E2 

-  5.411856  x  103E3  +  6.  17681 E4  x  105 

-  1.810118  x  107E5  -  6.2057  x  108E6 

+  4.406075  x  1010  E7  -  6. 58746  x  101 1 E8 

dF2/dE 


1 


E  x 
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Ej  -  0.0245 

2.49995  x  10'3  +  9.37472  x  10* *E 
-  4. 62461  E2  -  4.452203  x  10lE3 
+  3.751364  x  102E4 
dF2/dE 


FIG.  4  Continued 
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JSTAR  is  the  number  of  a  cell  just  ahead  of  the  shock.  Cell  No. 
160  lies  at  radius  4.0  cm.  M(J)  is  mass  of  J1*1  cell  in  grams, 
WH(J)  =  waste  heat  per  unit  mass  in  cell,  W{J)  -  total  waste 
heat  in  all  cells  up  to  and  including  J.  S ( J >  =  increase  in  entropy 
of  I1*1  cell  beyond  initial  value  SQ. 


FIG.  5  FLOW  CHART  FOR  CALCULATION  OF  WASTE  HEAT 
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V  NUMERICAL  RESULTS 


Values  of  S  -  S  ,  ,  and  W  are  given  in  Tables  V,  VI,  and 

VII  for  simple  water,  NOL  water,  and  simple  glycerine.  "Simple" 

means  with  constant  b  and  C  as  in  Section  II;  "NOL  water"  means 

v 

that  the  equation  of  state  of  Fig.  4  has  been  used.  J  is  the  number 
assigned  to  each  cell,  which  consists  of  a  spherical  shell  0.  02  cm  thick. 
The  initial  inner  radius  of  the  J  cell  is 

r j  —  1.6  +  0.  02  (J  -  41)cm  . 

The  mass  contained  in  the  J**1  cell  is 

4  (r  +  l)3  -  r j3  4r[rj‘  x  °*02 

Mj =  Tn  *470018 —  -  — rroorg—  grams 

The  mass  of  water  contained  within  and  including  the  J1*1  cell  is 

MJ  =  57TOt{[1*6  +  0-02(J-40)]3-(1.6)3J  (31) 

The  outer  boundary  surface  of  the  liquid  considered  in  this  example 
is  at  r  =  3.85  cm,  which  corresponds  approximately  to  the  inner 
boundary  of  the  153  cell.  Then  the  fraction  of  total  water  mass  con¬ 
tained  within  and  including  the  J**1  cell  is 


MJ  _  [  1.  6  +  0.  02  (J -40)]3  -  (1.6)3 
Mtotal  53 
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Table  V 
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WASTE  HEAT  AND  ENTROPY  FOR  SIMPLE  WATER 


J 

Entropy 

Change 

(Mbar  cc/g  deg) 

Specific 
Waste  Heat 
(Mbar  cc/g) 

Waste  Heat 

Inside  and  Including 
Jth  cell 
(Mbar  cc) 

41 

0.  000026 

0.  007620 

0.  004956 

51 

0.  000027 

0.  007914 

0.  071979 

61 

0.  000020 

0. 005861 

0.  133249 

71 

0.  000015 

0.  004396 

0.  189216 

81 

0.  000011 

0.  003223 

0.239713 

91 

0. 000009 

0.  002637 

0.286267 

101 

0. 000007 

0.  002051 

0.  328033 

111 

0.  000005 

0.  001465 

0. 366175 

121 

0. 000004 

0.  001172 

0.400461 

131 

0.  000004 

0. 001172 

0.432881 

141 

0.  000003 

0.  000879 

0.460213 

151 

0. 000002 

0. 000586 

0.487570 

153 

0. 000002 

0. 000586 

0.491905 

Table  VI 


WASTE  HEAT  AND  ENTROPY  FOR  NOL  WATER 


J 

Entropy 

Change 

(Mbar  cc/g  deg) 

Specific 
Waste  Heat 
(Mbar  cc/g) 

Waste  Heat 

Inside  and  Including 
Jth  Cell 
(Mbar  cc) 

41 

0.  000021 

0.  008026 

0.005219 

51 

0.  000021 

0.  007978 

0. 073345 

61 

0.  000016 

0. 005556 

0.  133261 

71 

0. 000012 

0.  004069 

0.  186032 

81 

0. 000009 

0.  003037 

0.232691 

91 

0. 000007 

0.  002317 

0.274368 

101 

0.  000006 

0.  001774 

0.  311459 

111 

0.  000004 

0.  001379 

0.  344459 

121 

0.  000004 

0.  001085 

0.373955 

131 

0.  000003 

0.  000866 

0.400460 

isi 

0. 000002 

0.  000697 

0.  424317 

in I 

0.  000002 

0.  000566 

0.  445878 

m 

0. 000002 

0.  000544 

0.449941 
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Table  VII 


WASTE  HEAT  AND  ENTROPY  FOR  GLYCERINE 


J 

Entropy- 

Change 

(Mbar  cc /g  deg) 

Specific 
Waste  Heat 
(Mbar  cc/g) 

Waste  Heat 

Inside  and  Including 
Jth  Cell 
(Mbar  cc) 

41 

0. 000014 

0.  005772 

0.  004730 

51 

0.  000016 

0.  006533 

0.  076166 

61 

0. 000011 

0.  004315 

0.  136227 

71 

0. 000009 

0.  003050 

0. 186773 

81 

0. 000007 

0.  002212 

0.23017C 

91 

0. 000005 

0.  001651 

0.267907 

101 

0.  000004 

0.  001257 

0.  301107 

111 

0.  000003 

0. 000972 

0.  330488 

121 

0.  000002 

0. 000763 

0.  356636 

131 

0.  000002 

0. 000604 

0. 379993 

131 

0.  000002 

0.  000484 

0.400910 

m 

0.  000001 

0.  000392 

0.419732 

0. 000001 

0. 000375 

0.  423265 

Table  VIII  gives  some  values  of  this  ratio,  along  with  values 
of  the  fraction  of  total  waste  heat  which  is  deposited  within  and  in¬ 
cluding  the  J**1  cell. 


Table  VIII 

DISTRIBUTION  IN  J  OF  LIQUID  MASS  AND  WASTE  HEAT 


J 

i 

Mt/M,  .  . 

J  total 

Fraction  of  Total  Waste  Heat 

Simple  Water 

NOL  Water 

Simple  Glycerine 

60 

.  09 

.26 

.29 

.  31 

80 

.  18 

.48 

.51 

.54 

100 

.  335 

.  66 

.69 

.71 

120 

.  54 

.  81 

.83 

.84 

140 

.  71 

.93 

.94 

.95 

152-1/2 

1. 00 

o 

o 

• 

1.  00 

1.  00 
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It  is  evident  that  in  all  three  liquids  the  major  part  of  the  energy 
dissipation  occurs  in  less  than  one -third  of  the  mass.  This  is  the  kind 
of  effect  one  would  anticipate,  since  the  energy  dissipated  in  shock 
compression  decreases  rapidly  with  shock  amplitude.  The  total  waste 
heats  in  the  three  cases  are: 

Simple  water  .493  Mbar  cc 

NOL  water  .449  Mbar  cc 

Simple  glycerine  .422  Mbar  cc 

The  total  energy  of  the  explosive  is  1.401  Mbar  cc,  so  these  numbers 

represent  35%,  32%,  and  30%  of  the  explosive  energy,  respectively. 

The  insensitivity  of  the  waste  heat  to  the  equation  of  state  of  the  fill  is 
particularly  impressive;  one  must  conclude  that  only  the  grossest 
features  of  the  equation  of  state  enter  into  wave  propagation  in  spherical 
geometry. 

Throughout  this  calculation  the  question  of  phase  change  has 
been  ignored.  As  long  as  the  pressure  in  the  material  remains  above 
a  kilobar  or  so,  this  is  probably  a  satisfactory  procedure.  However, 
specific  entropy  increase  in  cells  near  the  explosive  is  great  enough 
that  the  liquid  may  evaporate  when  pressure  is  reduced  to  atmospheric, 
depending  upon  its  specific  and  latent  heats.  The  fraction  of  material 
so  affected  is  estimated  by  locating  the  cell,  Jg  ,  at  which  the  specific 
waste  heat,  W^(J) ,  is  equal  to  hg  ,  the  heat  required  to  raise  unit 
mass  to  the  boiling  temperature  from  20°C.  This  heat  is  entered  in 
column  3  of  Table  IX. 

The  mass  of  liquid  contained  within  and  including  Jg  is  denoted 
Mg  and  is  calculated  from  Fq.  (31).  The  heat  required  to  boil  Mg 
at  atmospheric  pressure  is  (hg  +  Ly  )Mfi  ,  where  Ly  is  the  latent  heat 
of  vaporization.  The  waste  heat  within  and  including  cell  Jg  is  denoted 
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by  Wg  .  The  total  mass  of  liquid  which  will  evaporate  when  the  liquid 
expands  to  atmospheric  pressure  is  My  =  (Wg  -  h^Mg  )/Lv  .  Except 
for  decomposition  this  material  will  recondense  after  cooling  by  the 
atmosphere.  It  will  not,  however,  participate  in  the  liquid  expansion 
process  described  elsewhere  in  this  report.  It  will,  in  fact,  almost 
certainly  be  augmented  as  a  result  of  heating  through  direct  contact 
with  detonation  gases.  Some  values  of  My/M^^^  are  entered  in 
Table  IX.  These  have  been  calculated  on  the  assumption  that  the  distri¬ 
bution  of  waste  heat  for  each  liquid  listed  is  the  same  as  that  calculated 
for  NOL  water.  Table  VI.  Distributions  for  each  liquid  listed  in 
Tables  I  and  II  can  readily  be  calculated  according  to  the  procedures 
described  in  Section  IV,  but  the  differences  between  the  estimate  of 
Table  IX  and  a  more  precise  calculation  are  expected  to  be  small. 

Until  better  equations  of  state  are  constructed,  such  differences  are 
not  likely  to  be  very  significant.  In  order  to  illustrate  the  magnitude  of 
error  associated  with  the  estimating  procedure,  a  flow  calculation 
according  to  Section  IV  has  been  made  for  ether,  which  is  likely  to  ex¬ 
hibit  the  largest  deviation.  This  calculation  is  described  in  Appendix  A, 
and  the  results  are  entered  in  parentheses  in  Table  IX.  Since  ether  is 
the  most  compressible  of  the  liquids  in  Table  IX,  glycerine  is  the  least 
compressible,  and  water  is  intermediate  between  them,  these  three 
materials  should  bracket  the  behavior  of  most  liquids.  Thus  the 
values  in  the  last  column  of  Table  IX  can  be  taken  to  reasonably  repre¬ 
sent  the  effects  of  vaporization  for  various  liquids. 

The  last  column  of  Table  IX  suggests  that  experimental  com¬ 
parisons  of  distributions  in  ether  and  water,  for  example,  would 
certainly  indicate  whether  or  not  evaporation  has  a  significant  effect 
on  mass  distribution. 
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VI  CONCLUSIONS 


The  calculations  presented  here  show  that  the  equations  of  state 
of  liquids  developed  in  this  work  suffice  for  calculating  the  magnitude 
of  energy  absorbed  in  heating  the  liquid  fill  of  a  spherical  burster  to 
within  about  10%  of  its  correct  value.  For  the  particular  dimensions 
considered  here,  this  value  is  about  a  third  of  the  explosive  energy. 

The  explosive  is  similar  to  C-3and  the  ratio  of  liquid  mass  to  explosive 
mass  is  about  10.  Lower  explosives  would  produce  lea  s  heating  and  different 
explosive  configurations  couldbe  expected  to  alter  the  heating  effect. 

The  Q-code  used  to  perform  the  waste  heat  calculations  appears 
to  be  completely  satisfactory  for  the  purpose;  care  must  be  taken  to 
insure  that  cell  size  is  small  enough  to  yield  correct  values  of  the 
entropy  increase  in  shock. 

If  more  precise  values  of  the  waste  heat  are  required,  improved 
equations  of  state  will  be  needed.  Substantial  improvement  over  the 
present  equations  is  possible,  but  careful  thought  will  be  required  to 
make  optimum  use  of  available  thermodynamic  data.  Additional 
Hugoniot  points  would  be  of  considerable  help  in  constructing  an  im¬ 
proved  equation. 

If  serious  decomposition  of  fill  materials  is  expected  to  occur 
in  the  entropy  ranges  achieved  in  these  calculations,  the  calculations 
should  be  repeated  for  the  particular  materials  of  interest,  using  the 
best  equations  of  state  that  can  be  constructed.  This  is  because  specific 
entropy  (and  hence  temperature)  is  more  sensitive  to  thermodynamic 
properties  than  is  waste  heat. 

Although  calculations  have  not  been  made  specifically  for  pow¬ 
dered  solids,  it  is  anticipated  on  theoretical  grounds  that  waste  heats 
will  be  as  great  as  or  greater  than  those  calculated  for  liquids  because 
of  the  porosity  in  the  initial  fill.  Calculations  for  a  powdered  solid 
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can  be  made  using  initial  bulk  density  and  equation  of  state  of  the 
solid  material. 

The  amount  of  material  vaporized  by  passage  of  the  shock  differs 
markedly  between  water  and  other  materials  such  as  ether  and  alcohol. 
This  suggests  that  comparative  experiments  with  different  materials 
may  provide  information  about  the  importance  of  evaporation  in  dis¬ 
semination  by  an  explosive  burster. 

The  results  on  evaporation  also  suggest  that  this  may  be  a  factor 
influencing  the  relative  performances  of  bursters  containing  high  and 
low  explosives. 
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APPENDIX  A 


WASTE  HEAT  IN  ETHER 


The  calculations  described  in  Section  IV  were  performed  for 
ether  after  Table  IX  had  been  constructed  according  to  the  approxi¬ 
mations  described  there.  In  order  to  avoid  calculating  negative 
pressures,  as  shown  in  Table  IV,  a  set  of  coefficients  A,  B,  C,  D 
different  from  those  in  Table  III  was  used.  A^  of  Eq.  (23)  was  set 
to  zero,  A  was  determined  from  the  adiabatic  sound  velocity,  and 
Aj  was  determined  by  least  squares.  The  following  constants  were 
obtained: 


b  =  0.  5270 

Cv=  2.  26  x  10”5 

Tq=  293°K 

V  =  1.405 

o 

A  =  0.01511 

B  =  0.4428 

C  =  0.  1079 

D  =  0 

The  calculation  yielded  the  values  of  waste  heat  and  entropy 
shown  in  Table  A-l.  The  total  waste  heat  is  0.4045  Mbarcc/g,  and 
the  mass  of  material  vaporized  is  95.  5  grams,  or  0.  603  of  the  total 
liquid  mass.  These  numbers  are  shown  in  parentheses  in  Table  IX 
for  comparison  with  the  estimates  made  earlier.  The  mass  vaporized 
is  somewhat  lower  than  the  estimate,  but  the  difference  is  not  particu¬ 
larly  significant  in  view  of  the  approximations  made. 
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Table  A-l 


WASTE  HEAT  AND  ENTROPY  FOR  ETHER 


Entropy 

Change 

(Mbar  cc/g  deg) 

Specific 
Waste  Heat 
(Mbar  cc/g) 

Waste  Heat 

Inside  and  Including 
Jth  Cell 
(Mbar  cc) 

0.  000023 

0.  011758 

0. 005452 

0.  000020 

0. 009340 

0. 065224 

0.  000015 

0.  006360 

0.  114439 

0.  000012 

0.  004877 

0.  158461 

0. 000010 

0.  003717 

0.  198930 

0.  000008 

0.  002867 

0.235426 

0.  000007 

0.  002287 

0.268880 

0. 000006 

0.  001839 

0.299810 

0. 000005 

0.  001484 

0.  328244 

0.  000004 

0.  001215 

0.354422 

0. 000003 

0.001009 

0.  378729 

0. 000003 

0.  000844 

0.401349 

0. 000003 

0.  000814 

0.405678 
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